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We discuss the evaluation of certain d dimensional angular integrals which arise in pertur- 
bative field theory calculations. We find that the angular integral with n denominators can 
, be computed in terms of a certain special function, the so-called _ff-function of several vari- 

CN . ables. We also present several illustrative examples of the general result and briefly consider 

S : some applications. 

I. INTRODUCTION 

' When computing higher order corrections in perturbative field theory, the following d dimen- 

pi • sional angular integrals are encountered in many situations 
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where p^, . . . ^pn are fixed vectors in d = 4 — 2e dimensional Minkowski space and dD.d~i{q) is 
the rotationally invariant angular measure in d dimensions for the massless vector q^. For a 
single denominator, i.e. n = 1, the integral in Eq. ([1]) is easy to evaluate, as it reduces to a single 
(trivial) integration in a properly chosen Lorentz frame. The case of two denominators, n = 2, 
is already quite a bit more cumbersome, and it seems that general (i.e. ji and j2 are symbolic) 
analytic expressions valid to all orders in e are only available in the literature for the massless 
case {p\ = P2 = 0), as first derived in Ref. [ll]. When one or both of the momenta p'^ and 
massive. Appendix C of Ref. provides a very useful compilation of known results. (See also 
and references therein.) However, these are first of all limited to specific values of ji and j2, 
specifically ji,j2 = —2, —1, . . . , 2. (Some result^ for different specific values of j's — in particular 
integers with |ji|,|j2| < 4 — are also known [4].) Furthermore, they are given as expansions in 
e, up to and including 0(e) terms for the case of a single massive momentum (e.g. pf 7^ and 
P2 = 0), while for the case when both momenta are massive (pf 7^ and P2 7^ 0), only the four 
dimensional result is given. (Clearly the integral in Eq. ([T]) is finite in four dimensions if all pf, 
i = l,...,n, are massive.) Work towards deriving the 0(e) terms for the angular integral with 
two denominators and two massive momenta was presented recently in Ref. [5]. As explained in 
Ref. Q], the most difficult of these two denominator integrals were computed by relating them to 
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the imaginary parts of certain box integrals, which could be evaluated by Feynman parameters. 

However, in certain cases, results going beyond those found in Refs. 1 
example, when integrating [Q\ the so-called iterated singly-unresolved approximate cross section of 
the NNLO subtraction scheme of Refs. 7-[ll|, one requires a general (i.e. symbolic ji and 72), all- 



order (in e) expression for the two denominator angular integral with one massive momentum. In 
the same computation, one also encounters angular integrals with three denominators and general 
exponents. To the best of our knowledge, there is no systematic discussion of such d dimensional 
angular integrals with more than two denominators in the published literature. 



In this paper, we use the method of Mellin-Barnes representations (see and references 
therein) to evaluate the integral in Eq. ([T|) for arbitrary n, and massless or massive momenta 
(i = 1, . . . , n). The exponents ji {i = 1, . . . ,n) are also kept symbolic, and we tacitly assume that 
they satisfy any constraints that are needed to make our manipulations meaningful. In particular, 
it will be seen that our final expression for the general angular integral (and indeed its derivation) 
cannot be applied naively for nonpositive integer exponents. Nevertheless, some of our specific 
results will be valid even for ji being a nonpositive integer. 

The analytical expression for the general angular integral with n denominators is computed 
in Sect.|TIl and is given in terms of the //-function of several variables. The //-function of sev- 
eral variables has been discussed in various forms by a number of authors in the literature, see 
e.g. Ref. jisl] and references therein. (See also the recent book [l^, which however deals mostly 
with the single variable case.) For convenience, we recall the definition of the //-function as used 
in the present paper in Appendix El Then, in Sect. lIIH we illustrate the general case by several 
specific examples. In particular, we rederive and extend all known results for n = 2 as special cases 
of the general expression, including a general, all-order (in e) formula for the case with a single 



massive momentum. We also discuss the three denominator angular integral arising in Ref. 
draw our conclusions in Sect. lIVI 



. We 
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II. ANGULAR INTEGRAL WITH n DENOMINATORS 



A. General result 

To begin, we note that the overall normalization of the and plays no essential role, since 
clearly 
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Hence, it is no loss of generality to choose the normalization of all vectors in whatever way is most 
convenient. In particular, to write the integral in Eq. ([1]) explicitly, one may choose a Lorentz frame 
where 

p^ = (l,Od_2,/3i), 



P2 = (l,0<i-3,/32SinX2^\/32COSX2^^), 
Ai /I n ■ (2) ■ (1) n (2) . (1) n (1 

Ps = (l)0d-4,P3smx3 smxa % /Ss cos Xs sm Xg S/^scosXa 



(3) 



n-l 



n-2 



k=l 



k=l 



n-l 



n-2 



= (1> Od-l-n, Pn X^n^ , Pn COS X^n H • • • , Pn COS X^f^ sin xi^^ , /3n COS x[ 

while reads 

= (1, ..'angles'.., cos ■!?„ ]^ sint?fc, cos "dn-i sini?fc, . . . , cos i?2 sini?i, cos , (4) 

k=l k=l 

and we have used the freedom to choose the normalization to fix each zeroth component to be one. 
In Eq. Q, the notation ..'angles'., stands for the d — 1 — n angular variables that may be trivially 
integrated in Eq. ([1]). The explicit expression for the measure dTtd-iil) reads 



df^d-i(g) = n d(cosi?fc) (sint?fc)-^'+i-2^dOrf_i_„(g) , 



(5) 
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and hence Eq. ([T]) leads to the integral 
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which we take as the definition of Oj 

forms to that of Ref. while Refs. 0,131 use a different normalization (the factor of / dil.d-i-n{Q) is 



m=l 

j„. Notice that for n = 2 the normalization of f^jij2 (^on- 



4 



missing). As defined above, is a function of the "^"2 independent angles Xi \ • • • > XrT ^\ 

and n velocities, /3i, . . . However, it will be more natural to adopt the dot products between 
the various in Eq. as the independent variables. We will set the following notation 

(7) 



Vkl 



pI 



k = I 



where the choice of normalization will become clear later. If the pf are all light-like or time-like 
(i.e. < /3i < 1), then we have Vki > 0. In the following, we will assume that all Vki are nonnegative. 
We can now state our main result: the function ilj^ is given by the following expression: 



^h,-MM;e) = 2^-^-^'7T'-'H[v;{cx,Ay,iP,B);Ls] 



(8) 



where H is the i?- function of A'^ = "^"^^^'^ variables [13(. In Eq. ([8]) above, v denotes the vector of 



N variables, V = {vn,Vl2, ■ ■ • , f in, ^"22, ^23, • • • , fn-ln, Wnn) 

of parameters 



while a. and (3 are the following vectors 



« = (Otv, Ji, . . . , jn, 1 - j - e) , /3 = (ji, . . . , j„, 2 - j - 2e) 
In Eqs. ([8]) and dl]), j is the sum of exponents, 



(9) 



k=l 



(10) 



Notice that the number of components of cx is (!i±lK!l±Hl^ while that of /3 is (n + 1). Finally, A 
and B are (^^±i^x and (n + 1) x A matrices of parameters, respectively. We have 
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i.e. B is zero, while the n x N dimensional matrix M has the following block form: 
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In Eqs. and dill), 1 

axa denotes the a x a dimensional unit matrix, while (0)^x6 denotes an 
a X b dimensional block of zeros. To give some examples, we spell out the A matrix explicitly for 
the cases n = 1, 2 and 3, when A is a (3 x 1), (6x3) and (10 x 6) dimensional matrix, respectively: 
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(13) 

For clarity, in Eq. (jl3p . we have indicated the block structure of the various A matrices explicitly. 
Finally, in Eq. ([8]) we have Ls = Lg^ x . . . x Lg^^, where Ls^ is an infinite contour in the complex 
Sfc-plane running from — icxD to +ioo, whose properties we discuss below Eq. (fTSj) . Here x indicates 
the Cartesian product of contours. 

We note in passing that the /^-function of several variables satisfies various contiguous relations 



15l4l7l|. i.e. algebraic relations between functions H[v; (a, A); {(3, B); Ls] with the vectors of pa- 



rameters a and /3 shifted by vectors of integers. These relations may be used to reduce i/- functions 
to a set of basis functions with parameters differing form the original values by integer shifts via 
the method of differential reduction [3]. (See also Ref. [3] for a short but clear introduction to the 
main ideas.) The differential reduction of iZ-functions is beyond the scope of this paper. Neverthe- 
less, since the parameter of dimensional regularization, e, appears in ofsee Eq. ([9])), we speculate 
that this reduction will naturally include dimensional shift identities 201] for angular integrals. 

Finally, using the defining Mellin-Barnes representation of the /f-function as recalled in Ap- 
pendix [Aj we find 

^ji,...,ju{{vki};e) 



z vr 



nLir(jfc)r(2-i-2e) 
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Note that the Mellin-Barnes integration variables are denoted z^;, with k = l,...,n and 
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I = k,. . . ,n, i.e. zn, zu, ... , zin, ^22, Z23, • • • , Zn-in, Znn and in Eq. ([Tl|), we have furthermore 
introduced the following notation: 

ran k n 

z = ^ ^ Zfc, , and Zk = '^ zik + '^Zki. (15) 

k=l l=k 1=1 l=k 

In words, z is the sum of all N Mellin-Barnes variables, while z^ is the sum of all variables that 
involve k as one of their indices, such that z^k itself is counted twice, i.e. z^ = zik + . . . ^fc-ifc + 
2zfcfc + Zfcfe+i + • • • + Zkn. In Eq. p^ . j is the sum of all exponents, see Eq. p^ . The contours of 
integration for the Zki are chosen in the standard way: the poles with a,T{. . . + Zki) dependence 
are to the left of the contour and poles with a,T{. . . — Zki) dependence are to the right of it. 



B. Computation 



We establish Eq. ()14p by direct computation, as follows. Consider Eq. ([T|). First, use Feynman 
parametrization to combine all n denominators: 

-3 
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XkPk 
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where again j is the sum of exponents as in Eq. (jlOp . and the are given in Eq. ([3]). By rotational 
invariance, we can choose a frame such that 



'^XkPk = {^,0(1-2, and = (1, ..'angles'. ., sin ■!?, cost?) . 



(17) 
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Then we have 
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,^kPl = 'i^^XkXlVkl 
k=l l=k+l k=l k=l l=k 



1-/3^ = '^XkXi{pk-pi) + Y,- 
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In the frame of Eq. (I17p . the integral in Eq. (I16p reduces to 
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X / dl^d-2 / d(cost?)(sin7?)-^^[l-/3({xfc,t;fez})cosi?]-^ 
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where /3{{xk,Vki}) is given (implicitly) in Eq. (jlSp . Hence, we can exchange all but one angular 
integration for an integration over a Feynman parameter. 

The angular integral which appears on the second line of Eq. (jl9p is just the one denominator 
massive integral 17 j, to be discussed in more detail in Sect.[TTlBl For now, we simply state that it 
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has a Mellin-Barnes representation of the form 

0, = 22-^-2%i-^— -— i 

' r(i)r(2 - j - 2e) 



d£o 

27ri 



r(-zo)r(i + 2zo)r(i-j-e-zo) 



1 



(20) 



however, we defer the derivation of this until Sect.[TTlBl Using Eq. ([20]) in Eq. ([HD, we obtain 
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I- I3'^{{xk,vu]) 



Next, we perform the integral over the Feynman parameters. The only nontrivial x dependence 
and this is given by (see Eq. (fT8]l above) 
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which explains our choice of normalization in Eq. d?]). We can factorize all x dependence in Eq. ()2ip 
by writing ( ^~^ {{^k,t>ki}) \ 



as a multidimensional (in fact (A^ — 1) dimensional; recall that N 



n(n+i) ^ ]y[giii]2-Barnes integral: 
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Substituting Eq. (p3]) into Eq. ([2T]) . we obtain 
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Setting zq — Efc=i Er=fc ■^'s' = -^"^ changing the variable of integration zq — >• Znn, we find 
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r(j + 2z)r(i-j-e-z) 
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where z is the sum of all N integration variables as in Eq. (jlSp . Collecting all factors of the various 
x's in Eq. (I25p . we obtain 
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where Zk is defined in Eq. (|15p . We can now perform the Feynman parameter integrals via 



-'^ k=i \k=i / r 



l^k=lPk 



(27) 



to obtain Eq. (jHl), as claimed: 
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In writing Eq. (p8]) . we used X]fc=i(ifc + -^fe) = J + 2z. This completes the calculation. 

Before moving on, some comments are in order. First, note that the derivation of Eq. (j28p 
implicitly assumes that the exponents, jk {k = l,...,n) are not zero or negative integers, and 
indeed Eq. (128p is clearly not applicable as it stands when any jk is a nonpositive integer. In such 
cases, when e.g. —jk' G N, we can attempt to analytically continue Eq. ([28]) to the required value 
of jk', say by setting jk' — >• jk' + S and performing the analytic continuation 5^0. The analytic 
continuation of Mellin-Barnes integrals has been automated in the Mathematica package MB.m 



Second, we call attention to the fact that that Eq. (j28p was obtained under the assumption that 
Vki > for all k = 1, . . . ,n and / = k, . . . ,n. However, sometimes it may happen that some Vki 
is identically zero, as e.g. when, say, momentum p^, in Eq. ([3]) is massless, implying Vk'k' = 0. In 
such situations we clearly cannot use Eq. ()28p as it stands. Nevertheless, the derivation is trivial 
to adapt to such cases, since when say Vk'i' is identically zero, the only change is that this term 
is missing form the sum in Eq. (|22p . Then, the corresponding Mellin-Barnes integration over Zk'i' 
is absent in Eq. ()23p . but the rest of the derivation goes through unchanged. The end result is 
that we must drop integrations corresponding to variables that are identically zero from the final 
expression, Eq. (I28p . Ultimately, this amounts to simply restricting all products (such as the one in 
the first bracket on the second line of Eq. ([28]) ) and sums (as in the definitions of Zk and z, Eq. (fT5|) ) 
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over Zf^i to values of k and I such that v^i 7^ 0. Needless to say, the i^-function representation of 
the integral must also be adapted to accommodate the fact that some integration variables are 
missing. 

III. EXAMPLES 

In this section we illustrate the use of the general result in Eqs. ([8]) and (JH]) with several 
examples. 

A. One denominator, massless 

We begin with the simplest example, the massless one denominator angular integral, i.e. n = 1 
and Pi = (hence /3i = 1). In this case, Eq. ([6]) reduces to 

nj{0;e) = j dnd~2 j d(cos??i)(sinT9i)-2^(l -cosT?i)~-''. (29) 

Since is massless, vn is identically zero, and the discussion at the end of Sect. Ill Bl applies. Thus, 
the Mellin-Barnes integral representation is zero dimensional and so clearly zi = z = 0. Then we 
find 

QM e) = 22-i-2%i-iXLzl^ . (30) 
' r(2-j-2e) ^ ' 

The result in Eq. ()30p is easy to verify by explicit computation: recalling that 

27r2 



and setting cosi?i — )• 2s — 1 in Eq. (j29j) . we obtain Eq. (j30j) immediately. 



B. One denominator, massive 

The next simplest example is the massive one denominator angular integral, i.e. n = 1, but 
p\ 7^ 0. In this case, Eq. ([6]) gives 

nj{vn;e) = J dnd-2 J d(cos??i)(sin??i)-2^(l-/3icos??i)-^ (32) 



Now vii is nonzero, and Eq. (|14p yields the following one dimensional Mellin-Barnes integral rep- 
resentation 

1 r-+ioo J 

"^'"'"^) = ^"''"-'" r(;)r(2-j-2.) /,„ ^r(-.„)r(; + 2.„)r(i-,-.-.„)(»„)- . 

(33) 
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We used that Eq. ()15p gives zi = 1z\\ and z = z\\. In terms of the i?- function, we have 

^■(m;e) = 22-^-2V-^/7[(t;n);(«,A);(/3,B);Ls], (34) 
where q, /3, A and B are given in Eqs. ([9|)- ()13p . with n = 1. ExpHcitly 



« = (0,j,l-j-e), P = {j,2-j-2e), 



(35) 



and 



-1 



S = [(0)2x1] 



(36) 



We may compute the integral in Eq. (|33p by using the doubhng relation for the gamma function, 

22X-1 



T{2x) 

to write Eq. (|33p in the following form 



-r(x)r(x + i) , 



(37) 



1 



T{j)T{2-j-2e) 



dzu 



(38) 
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r(-zn)r (I + r ( + zn r(i - j - 6 - zii)(4^ii)^" 



Then the Mellin-Barnes integral on the second line can be evaluated in terms of T functions and 
the 2F1 hypergeometric function. Indeed, we have 

r(c) f+'°° dz . 



2Fi{a,b, c,x) 



T{a)T{b)T{c-a)T{c-b) 



2m 



T{a + z)r{b + z)r{c -a-b- z)r{-z){l - xY , 



(see e.g. appendix D of Ref. [12]), and hence we find 



(39) 



nj{vii;e) = 2^-^'7r 



r(j)r(2-i-2e) 

r(i)r(m)r(3^-e)r(^-e 



(40) 



r(|-e) ' '\2' 2 '2 

Using Eq. (|37p . we can clean up the prefactor and obtain the final expression 

n,(.„;.) = 2-..-i|^.F.(|,^,^-,l-4„„) . (41) 

This result is also simple to verify by explicit calculation, since Eq. ()32p is straightforward to 
evaluate via the substitution cost?i — >• 2s — 1. We obtain 



n,{vu;e) = + P^r^F, [j, l - 6,2 - 2e, ^ 



(42) 
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and Eq. (jH]) is then reproduced using the quadratic hypergeometric identity (see e 

2n 



•g- 



2Fi{a,b,2b,z) = \Fi 



a a + 1 , 1 

,b+ -, , 

2' 2 ' 2' V 2 - z 



(43) 



and the relation vu = — The above considerations in fact estabUsh Eq. (jSSp independently of 
Eq. ()14p . thus the gap left in the derivation of the main result is closed. 



Before moving on, we note that although Eq. (jSSp was derived under the assumption that j is 
not zero or a negative integer, our final result, Eq. ()4ip . is in fact valid for negative integer j as 
well. 

In practical applications, one is often interested in the e-expansion of the final result, Eq. (j4ip . 
When j is an integer, it is straightforward to obtain such expansions starting from the equivalent 
form of the result, Eq. (j42p . Indeed, for j a negative integer, the power series representation of 
the hypergeometric function in Eq. (|42p terminates, and the e-expansion in this case is trivial. For 
positive integer j, the method of nested sums of Ref. or the integration method of Ref. [2^] may 



be employed. The nested sums method has been implemented in several publicly available packages 



such as nestedsums [25], XSummer 26| and HypExp {2^ . with the last of these implementing the 
integration method as well. Algorithms have also been developed for the expansion of (generalized) 



hypergeometric functions around half-integer values of the parameters 
package 
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32l |. and the HypExp 2 



32( 1 provides a public implementation of one particular method. Finally, we note that 
whenever j is not zero or a negative integer, the direct numerical integration of the Mellin-Barnes 
representation in Eq. ( 1331) provides a fast and reliable way to obtain numerical results. 

on 

By way of illustration, and for purposes of comparing with existing literature [3) y], we obtain 
the e-expansion of Eq. ([12]) for the specific values of j = —2, —1, 1 and 2, up to and including O(e^) 
terms. (Note that for j = 0, ^o{v; e) just reduces to the massless integral f^o(0; e), and we do not 



discuss this case further.) The results, obtained with the method of nested sums and XSummer 
are presented in Appendix IC 1[ 



C. Two denominators, massless 

Our next example is the massless two denominator angular integral, i.e. n = 2, with pf = P2 = 
(hence /3i = /32 = !)• Eq. Q reads in this case 



nj^k{vi2,0,0;e) = J dnd-3 y^d(cos^9i)(sin??i)-2^y"^d(cos^92)(sin??2)^ 



-l-2e 



(44) 



X (1 — cos ^i) ■' {I — COS x^2^ COS I?! — siu x'2' sin cos 'i?2 



(1) 
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Since both p'^ and P2 are massless, we have vu = V22 = 0. Then the discussion at the end of 
Sect. lIIBl apphes and we find that Eq. ()14p leads to the following one dimensional Mellin-Barnes 
integral representation 

1 



^>(t'i2,0,0;e) = 2 



r(j)r(A:)r(2 -j-k-2e) 



-+1CO J 

X / —^T{-zi2)T{j + zi2)T{k + zi2)r(l -j-k-e- Z12) {vi2Y'^ . 



(45) 
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We used that Eq. (|T5|) gives zi = Z2 = z = Z12. The ff- function representation of Eq. (|45|) reads 



^■,fc(^^i2,0,0;e) = 2 



2-}-fc-2e„l-e 



7r'-'H[{vi2);ia,A);if3,B);Ls], 



(46) 



where 



a = (0, j, A;, 1 - i - A; - e) , /3 = {j, k,2-j-k-2e), 



(47) 



and 



-1 



B = [(0)3x1] 



(48) 



Notice that A above may be obtained from the general expression for A{n = 2) in Eq. (jl3p by 
removing the first and third columns which would correspond to the variables vu and V22 which 
are identically zero, and then removing the first and third rows of the matrix so obtained, which 
contain only zeros. Correspondingly, the first and third components of cx (both zeros) are also 
removed as compared to the general formula for n = 2 in Eq. ([9]). 

The Mellin-Barnes integral in Eq. (I45p straightforwardly evaluates in terms of T functions and 
a 2F1 hypergeometric function, see Eq. (p9]) . and we find 



n,,,{vuA 0; e) = 2^-^-^-^-vri- X^^^ ^^^^^f^j ,^ £ 2^1 {j, k,l-e,l- vu) 



r(l-e)r(2-j-fc-2e) 



(49) 



Upon noting that 



V12 



{Pi -^2) _ 1 - COSX2 



(1) 



1 - Vl2 



1 , (1) 2 (1) 

1 + COS X2 COS X2 



(50) 



the result in Eq. (|49p is seen to coincide with Eq. (A. 11) of Ref. [l|. 

We remind the reader that Eq. ()45p was derived under the assumption that j and k are not zero 
or negative integers. Nevertheless, the final result in Eq. ()49p applies in such cases as well. 
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Finally, we mention that the expansion of Eq. ()49p in e for integer or half-integer j and k 
is straightforward, as discussed at the end of Sect. lIIIBl Here, by way of illustration, we present 
these expansions for j, k = —2, —1,1 and 2 in Appendix IC 21 (We do not consider cases where either 
exponent is zero, since these are not genuinely two denominator angular integrals.) Expansions 
of the appropriate hyper geometric functions were computed with the nested sums method and 
XSummer M. 



D. Two denominators, one mass 



Now consider the generalization of the previous example to the single mass case, i.e. when, say, 
/ but pI = (hence /3i / 1, but /Sg = 1). Then Eq. dSD gives 



^j,k{vi2,vu,0;€) = I dO,d-3 j d(cosi9i)(sin-i?i) j d(cosi92)(sin'!92)" 



-l-2e 



(51) 



X (1 — /3i cos -iJi) ■' (1 — cosx2^^ cosT?i — sinx2^^ sini?! cost?2) • 
Since ^22 still vanishes identically, the Mellin-Barnes representation of Eq. (|14p is only two dimen- 
sional 

1 /■+'°°dzndzi2, 



%,fc(^^i2,wii,0;e) = 2 



T(-zn)r(-zi2) 



r(i)r(A:)r(2 -j-k-2e) J_,^ (2^i)2 
X r(j + 2zu + A2)r(fe + zi2)r(l -J-k-e-zu- Z12) {vnY^^ (vuT'' 



(52) 



We used that in this case, Eq. ([T5|) evaluates as zi = 1z\\ + Z12, z^ = z\2 and z = z\\ + zi2- Written 
in terms of the if-function, Eq. (j52p is 



^j,kivi2,vii,0;e) = 2 



2-j-k-2e l-e 



7r'''H[{vn,vi2y,{oc,Ay,{f3,By,Ls] 



where 



and 



a = {0,0,j,k,l-j-k-e), 13 = (j, k,2-j-k-2e) 



(53) 



(54) 



-1 

-1 

2 1 

1 

-1 -1 



B = [(0)3x2] 



(55) 
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The A matrix above is obtained from the general expression for A(n = 2) in Eq. (jl3p by removing 
the third column corresponding to the variable V22 which is identically zero, and then removing 
the third row of the matrix obtained, which contains only zeros. Accordingly, the third component 
of a (again zero) is also dropped as compared to the general formula for n = 2 in Eq. ([9]). 

The two dimensional Mellin-Barnes integral in Eq. ()52p can be evaluated in terms of the Appell 
function of the first kind. We show this in Appendix [Bl and only quote the final result here. We 
find 

%,,(.12,t^ll,0;6) = 2^-.-^-2e,l-e ra-fc-^^)^ 

2^12 - 1 - VI - 4fii 2vi2 - 1 + VI - 4^11 



X Fi ( j,l - k - e,l - k - e,2 - k - 2e 



2vi2 



2v 



12 



(56) 



Let us make several comments. First, as in the previous two examples, the final expression 
in Eq. (156^ is valid even for j or k zero or a negative integer, even though the Mellin-Barnes 
representation in Eq. ([52]) . as it stands, does not apply to these cases. 

Second, the Appell function of the first kind is precisely the type of generalized hyper geometric 



function w 
ofRefs. 231 



rose expansion around integer or half-integer parameters can be solved with the methods 



281432]. In Appendix IC 3( we present the expansion of Eq. (f56]l in e up to and including 
finite terms, for j,k = —2,-1,1 and 2. (Again, we only deal with cases which genuinely involve 
two denominator integrals.) We used the method of nested sums and XSummer [26] to compute the 
expansions of the appropriate Appell Fi functions. 

Finally, let us briefly discuss the application of the result in Eq. (|56p to the computation of 



As explained in 



certain integrated counterterms in the NNLO subtraction scheme of Refs. 
6|], when computing the so-called integrated iterated singly-unresolved approximate cross section, 
several integrals must be evaluated which involve the the function J'^^'^\Y, l3;e,yo,d'Q) in their 
integrands, where (see (C.25) of Ref. [s]) 



^(i-)(y,/3;e,yo,do) = -4y; 



r2(i 



-nn{cosxiY,/3) 



ryo 

,/3,l)/ 
Jo 



dyy' 



-l-2e 



(1 - ?/)"o 



(57) 



27rr(l - 2e) 

and Y as well as f3 depend on further integration variables. Because of this, we require an all- 
order (in e) result for J^(i"^)(y, /3; e, yo, 4)- In Eq. ([57]), Oii(cos x(^, /3), /3, 1) is a special case of the 
general function $7jfc(cos /32); defined in (C.19) of Ref. [g] as follows 



%fc(cosx,/3i,/32) 



d(cost?)(sin??) 



-2e 



d(cos 99) (sin ip)' 



-l-2e 



(58) 



X (1 — /3i cost?) ■'[1 — /32 (sin X sin cos (/9 -|- cos x cost?)]' 
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For = 1, which is the case relevant in Eq. (j57p . this is clearly just proportional to the one 
mass, two denominator angular integral. Hence the results of this section can be used to evaluate 
j(i-)(y,/3;e,yo,rfo) analytically. 



E. Two denominators, two masses 

Next, consider the general massive two denominator angular integral, when both pf ^ and 
P2 / (hence /3i 7^ 1 and /32 / 1). In this case, Eq. ([6]) gives explicitly 



j dCld-s j d(cos7?i)(sini?i)-2^y" d(cos7?2)(sini?2)' 



l-2e 



(59) 



X (1 — /3i COS!?!) ■'(I — /32 C0SX2 "^osT?! — /32 sinx2 sill"!?! COSt?2) ■ 
Now Eq. (|14p leads to the following three dimensional Mellin-Barnes integral representation 



^•,fc(^^i2,^^ii,'y22;e) = 2 



2-j-k-2e^-€ 
J Yj- 



1 



r(j)r(A;)r(2-i-A:-2e) 

dziidzi2dZ22T.. x-p. . 

— j^^^ — r(-zii)r(-zi2)r(-2;22) 

X r(j + 2zn + zi2)T{k + zu + 2^22) 



(60) 



X r(l - j - A: - e - zii - Z12 - ^22) (^^ii)"" (^^12)"^^ (^^22)^^^ • 

We used that Eq. (fT5]) gives zi = 2zii + Z12, Z2 = z\2 + 2z22 and z = z\\ ^ z\2 + ^22- Since all 
variables {y\\^ v\2 and V22) are now different from zero, the iJ-function representation of Eq. (j60]) . 

^jM^v^2.vxx.V22\ e) = 2^-^ -'=-2^1-^17 [(^;n, ^^12, t^22); («, A); (/3, S); Ls] , 



is simply the general expression in Eqs. (i9|l- (113p for n = 2, i.e. we have 

a = (0,0,0,j,/c,l -i-fc-e), /3 = (j, /c, 2 - j - 2e) , 

and 

'"-10 
0-10 
0-1 
2 1 
1 2 
-1 -1 -1 



(61) 



(62) 



B = [(0)3x3] . 



(63) 
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In this case, we are no longer able to evaluate the Mellin-Barnes integrals in Eq. (jGOp in terms 
of functions more familiar than the //-function of several variables. 

Nevertheless, the Mellin-Barnes representation of Eq. (160p is still a very useful starting point 
for computing the e expansion of ^(f 12, wii, ^22; e)- Let us briefly review the main steps involved. 



1. In general, the contours of integration in Eq. (j60p are not necessarily straight lines, and their 
standard definition is such that the poles with a r(. . . + z^i) dependence are to the left of 
the contour for the zm integration, while poles with a r(. . . — z^i) dependence are to the 



33| realized straight line contours parallel 



right of it. However, as a key observation, Ref. 
to the imaginary axis in an algorithmic way. The basic idea is that starting with a curved 
contour that fulfills the condition on the poles, one may deform it into a straight line, taking 
into account the residua of the crossed poles according to Cauchy's theorem. This procedure 
lends itself to implementation in computer codes for the evaluation and manipulation of 



Mellin-Barnes integrals, such as the MB.m package of Ref. 



2. Upon deformation of the curved contours, all potential singularities in e are extracted so 
that it is safe to expand in e around zero before performing the complex integrations. In this 
way, the Mellin-Barnes representations of the required coefficients of the Laurent expansion 
of the original integral are obtained. 

3. In the next step, we convert the complex contour integrations into sums over residua using 
Cauchy's theorem. 

4. Finally, we evaluate the sums. 

For the specific case of %,fc(fi2 5 I'll) ^'22; e) with integer j and k (we considered the cases j,k = 
—2, —1, 1 and 2 as before), this procedure leads to a representation of the 0{e^) coefficient which 
involves at most single sums. These are all straightforward to compute and we present the results 
in Appendix IC 41 

However, starting from the linear term in e, we are lead to representations of the coefficients 
involving up to triple sums, which are difficult to compute, and we made no severe effort to cal- 
culate them. In fact, this corroborates the findings of Ref. jsl nicely, where a completely different 
method leads to a one dimensional real integral representation of the linear term in the e expansion 



of nj^fc(vi2, uii, U22; e) which ". . . invo 
tion analytically. . . " and hence Ref. 



ves three square roots so it is difficult to evaluate the Integra- 



■ fl 



does ". . . not have an analytical answer." Nevertheless, we 
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remind the reader that the direct numerical integration of the MeUin-Barnes representation pro- 
vides a convenient and efficient way of obtaining numerical results for the higher order coefficients. 

F. Three denominators, massless 

As a final example, we present explicitly the massless angular integral with three denominators, 
i.e. n = 3 and pf = (hence f3i = 1) i = 1, 2, 3. Eq. ([6]) reads 



d(cos t?2)(sin 'd2) 



^2, ' 



^j,k,ii'"i2,vi3,V23;e) = j dild-A j d(cost?i)(sinT9i) "^^ J < 

/I 
d(cosi93)(sin??3)^^^^'^ (1 — cos'i?i)~-^(l — cosx2^^ cosi?! — sinx2^^ sini?! cos??2 

/I (1) a (2) ■ (1) ■ a Q . (2) . (1) . a . a q \~l 

X (1 — COS Xs COS Vi — COS X3 sm X3 sm vi cos V2 — sm X3 sm X3 sm vi sm V2 cos V3) 

(64) 

Since va = all for i = 1, 2, 3, the Mellin-Barnes representation of Eq. (I14p collapses to only three 
integrations 

nj,kA^i2,vi3,V23;e) = 22-J-'^-'-2^7ri-^ 



r(i)r(A:)r(/)r(2-j-A:-/-2e) 

/"+'~ dzi2d2i3dZ23T.. N-p, X 

^ / — — r(-zi2)r(-zi3)r(-z23) 

X r(j + Zi2 + Zi3)T{k + Z12 + Z23)r(/ + Zl3 + ^23) 
Xr{l-j-k-l-e-Zi2-Zi3- Z23) M''' ivi3r'' {V23r'' ■ 

Eq. (USD gives zi = Zi2 + Z13, Z2 = 212 + -223, -^3 = -^13 + -^23 and z = Z12 + Z13 + Z23 in this 



case, which we used when writing Eq. (|65p . In terms of the //-function, Eq. (|65p has the following 
representation 

%,fc,«(^^i2, ^^13, ^^23; e) = 2^-^-''-'-^'7r^-'H[{vu,vi3,V23y, (a, A); {f3, B), Ls] , (66) 

where 



a = {0,0,0,j,k,l,l-j-k-l-e), = {j,k,l,2 - j - k - I - 2e) , (67) 
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and 



-1 











-1 











-1 


1 


1 





1 





1 





1 


1 


1— ' 


-1 


-1 



B 



m 



4x3j 



(68) 



We may obtain the A matrix above from the general expression for A{n = 3) in Eq. (|13|) by 
removing the first, fourth and sixth columns corresponding to the variables Vn, V22 and ^33 which 
are identically zero, and then removing all rows of the resulting matrix, which contain only zeros 
(i.e. rows one, four and six). The corresponding components of a (i.e. the first, fourth and sixth, 
all zeros) are also dropped as compared to the general formula for n = 3 in Eq. ([9]). 

As in the previous example, we are again unable to evaluate the Mellin-Barnes integrals in 
Eq. (j65p in terms of functions other than the ff-function of several variables. 

However, as discussed in Sect. lIIIEl Eq. ()65p provides a useful starting point for obtaining the 
e expansion of f^j,A;,/(t'i2i ^i3i ^23; e)- In this particular case, for j, k and / integers (in fact we 
consider only j,k,l = 1 and 2), the procedure outlined below Eq. (j60p leads to zero dimensional 
Mellin-Barnes integral representations for both the 0(e^^) and 0{e^) coefficients in Step 2. Hence, 
there are no integrals or sums to compute at all. The expansions obtained, up to and including 
0{€^) terms, are presented in Appendix IC 51 

The situation with higher order expansion coefficients is very similar to the previous example 
of Sect. iniEl Again, starting from the linear term in e, we find representations involving triple 
sums which are hard to compute. In passing, we note that one may also attempt to evaluate the 
Mellin-Barnes representations of given expansion coefficients by means other than converting them 
into sums e.g. with methods along the lines of Ref. [3^, where the authors compute a difficult three 
dimensional Mellin-Barnes integral in terms of Goncharov polylogarithms. If we are satisfied with 
numerical results for higher order expansion coefficients, then the direct numerical integration the 
Mellin-Barnes representation proves convenient. 

We finish by noting that the knowledge of Eq. ()65p was necessary to compute certain iterated 
singly-unresolved integrals in Ref. [s]. Specifically, when computing the so-called integrated soft- 
double soft counterterm, we encounter the massless angular integral with three denominators in 
intermediate stages of the calculation. In particular, in order to be able to write the Mellin-Barnes 
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representation of the I^.Hjj^ master integral of Eq. (E.52) of Ref. Q], we required a Mellin-Barnes 
representation for the angular integral 

f d(cost?)(sini?)-2^ / d(cos (/?)(sin (1 - cos 

7-1 J-l ^gg^ 

X (1 — COS X2^^ cos 1? — sin ^2^'* sin cos (p)~^^ (1 — cos Xs^^ cos {} + sin xi'^ sin cos i/?)^^^ , 

where zi and 22 are integration variables of further Mellin-Barnes integrals. The integral in 

Eq. (j69|) is clearly proportional to ^^1,21,22(^12, "WiSi ^^23; e), in the special case where sinxs = and 
(2) 

cos X3 = — 1 in Eq. ([3]) . (This constraint means that only two variables out of f 12 , ^13 and ^23 are 
independent.) The results of this section thus provide the necessary Mellin-Barnes representation 
of Eq. ([69|) , and hence are needed to compute the j^. master integral. 

IV. CONCLUSIONS 

In this paper, we have evaluated some d dimensional angular integrals which arise in perturba- 
tive field theory calculations. We used the method of Mellin-Barnes representations to compute 
the general angular integral with n denominators, massive or massless momenta and (essentially) 
arbitrary powers of the denominators in terms of a certain special function, the so-called /f-function 
of several variables. We pointed out that the existence of various contiguous relations for the H- 
function provides the opportunity to apply the method of differential reduction to angular integrals 
in d dimensions. It would be very interesting to expand the present results in this direction. 

We illustrated the use of our general result with several examples of angular integrals with 
up to three denominators. We showed that some of these integrals can be computed in terms of 
(generalized) hyper geometric functions. In particular, the single denominator massless integral is 
fully expressed by F functions, while the massive integral involves the 2F1 hypergeometric function. 
For the massless two denominator integral we recover the known result of Ref. [l| which again 
involves a 2-^1 hypergeometric function. However, our derivation is much more straightforward 
than the original computation. When precisely one of the momenta is massive, we obtain a new 
all-order (in e) analytical expression for the two denominator angular integral which involves the 
Appell function of the first kind, Fi. 

In some applications, one is interested in the expansion of the angular integrals in the parameter 
of dimensional regularization, e. We discussed briefly how such expansions can be obtained starting 
from the corresponding Mellin-Barnes representations. By way of illustration, we have explicitly 
presented such expansions for all examples discussed, for a few specific values of exponents. 
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The results of this paper have already found applications in computing certain phase space 
integrals which appear when integrating NNLO subtraction terms. In fact, all specific angular 
integrals that are encountered in the integration of the so-called singly-unresolved and iterated 
singly-unresolved subtraction terms of Refs. ^, ^ were discussed in this paper explicitly. We 
expect that our present results will also prove valuable when computing the so-called integrated 
doubly- unresolved subtraction terms. 



Acknowledgments 

It is a pleasure to thank M. Y. Kalmykov, S.-O. Moch and Z. Trocsanyi for useful discussions and 
comments on the manuscript and E. Laenen for providing copies of Ref. and some unpublished 
notes of W. L. van Neerven. This work was supported in part by the Deutsche Forschungsgemein- 
schaft in SFB/TR 9 and. 



Appendix A: The i7-function of several variables 

In this Appendix, we recall the particular definition of the ff-function of several variables which 
we use in Sect. Ill Al This function has been discussed in various forms by several authors in the 
literature, here we adopt (essentially) the definition of Ref. ^]. In the most general case, the 
i?- function of variables is defined as follows: 

H[x,{a,A),{(3,B);Ls] = {2TT[)-^ f e(s)a;^ds, (Al) 

J Lo 



= 7- — f • (A2) 



where 



Here s = (si, . . . , sat), x = (xi, . . . ,xn), ol = (ai, • • • , CKm) and /3 = (/3i, . . . , /3„) denote vectors of 
complex numbers; while 

^ = {(l'j,k)mxN and B = {bj^k)nxN (A3) 

are matrices of real numbers. Also 

N N 



X^ 



YlixkY"; ds = JJdsfc; = x . . . x L^^ , (A4) 



k=l k=l 
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where Lg^ is an infinite contour in the complex s^-plane running from — ioo to +ioo such that 0(s) 
has no singularities for s G Ls- 

The ff-function of Eq. ()A1|1 generalizes nearly all known special functions of N variables, 



e.g. Lauricella functions F^^'^ and F^^; the G- function of N variables; the spe- 

cial i7-function of N variables, etc. For the specific cases of = 1 and 2, it essentially reduces to 
the known Fox's ff-function of one variable and the ff-function of two variables defined by various 
authors scattered in the literature. The definition given in Eq. (lAll) is different form the i7- function 
considered by Ref. [13] only in the replacement of £C ^ by a;^. We have made this replacement for 
convenience in our applications. 

Appendix B: A Mellin Barnes Integral 

In this Appendix, we evaluate the following two dimensional Mellin-Barnes integral analytically 

/+ioo J J 
-^r(-zi)r(-z2)r(a + 2zi + z^nt + Z2)r(c - zi - ^2) x^^ y^^ , (Bl) 

which we encounter in Sect. lIIIDl Throughout this Appendix, we assume tacitly that all parameters 
and integration variables lie in a strip of the complex plane such that each integral we write 
converges. 

We begin by writing the product of the second and fourth gamma functions above as a one 
dimensional real integral: 

/ = r(6) [' dss-'-'Hl-s)'-'+^'T{-zi)r{a + 2zi + Z2)T{c-zi-Z2)x^^y^^ . (B2) 

J-ioo (2vri)^ Jo 

Now, it is somewhat easier to follow the manipulations below if we make the change of variables 
Z2 — > —a — 2zi — Z2. 

'■+'°°dzidz2 /"^ j„ „a-l+2zi+Z2/i „\fe-a-l-2zi-Z2 



-y \ a+2zi+Z2 



(B3) 



X r(-zi)r(-z2)r(a + c + zi + 22) - 

\y 

Next, we rearrange some factors and write / in the following form 

/ = y"«r(6) rdss''-i(l - s)"+''+2'=-i 
Jo 

/■+^°°dzidZ2^, , fxs^y \s(l- S)Y\, 



\21 -~a.-c-zi-Z2 

s) \ 



(B4) 
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The Mellin-Barnes integrals are now easy to perform and we find 



I = y-''T{b)T{a + c) / dss^-^l-s 
Jo 



.a+b+2c-l 



s(l — s) xs 



+ 



Factoring the quadratic expression in the square brackets, 



s(l — s) xs 



+ 



22/ 



^ _ 2y-l- Vl-4x \ / _ 2j/-l + Vl-4x ^ 



2y 



we obtain 



I = y-^r{b)r{a + c) [ dss^-^l-s) 
Jo 



a+b+2c-l 



27/ - 1 - VT^^ 
X 1 s 



^ 2y-l + TT^l^ 
1 ; S 



(B5) 



(B6) 



(B7) 



2y J \ 2y 

The final integral can be performed in terms of the Appell function of the first kind (see e.g. 22l |) 
and we find 

„„ r(a)r(b)r(a + c)r(a + b + 2c) 
~^ r(2a + 6 + 2c) 

X Fi ^a, a + c, a + c,2a + b + 2c, 

which is our final result. 



2y - 1 - VI - 4x 2y - 1 + Vl - 4x 



2y 



2y 



(B8) 



Appendix C: Expansions 



1. One denominator, one mass 



In this section, we present the e-expansion of $7j(v;e) for the specific values of j = —2,-1,1 
and 2. In order not to clutter the following expressions with irrelevant constants like ln(47r) and 
7£;, we extract a factor of J dild-s and present the e-expansion of the function Ij, where 



w,. ^ _o-l+2. . r(l-26) 



%(w;e) , 



(CI) 
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which has the further advantage that it may be directly compared with the appropriate expressions 
of Refs. y,y. We find 



I-2(v;e) 



[4 4 /26 32 \ /160 208 \ n , 

2tt\ v + i v]e+ { V + 0(e^ 

[3 3 V9 9y V 27 27 J ^ 



l + 2e + 4e2 + 0(e3) 



(C2) 



(C3) 



hiv;e) 



IT 



VI -^v\ VI- Vi-4vy 2 



1 - y/l-Av 



i + 7r^4^ 



1 

+ 6 



In 



3 j 1 + yi - 4t; 
1 - VI - 4t; 



1 + vr^4^ 



1 - VI -4u 



+ 24 Lis 



1 - VI - 4t; 
2V1 - 4t; 



1 + vr^4^ 



1 + VI - 4vy VI - VI -4t; 

1 - VI -4v 



12 Lis 



1 + VT^^ 



12 C: 



/2 ft;; e) 



vr 



4WV1 - 4-u 



2Vr^+21n(l±^^U + 

' 1 - vr^4^ 



In 



1 + vr^4^. 

e' + 0(63)|, (C4) 
2 / 1 + VT^^^ 



1 - VI -4f 



+ 4 Lis 



2V1 - 

1 + VI - 4t; 



(C5) 



2. Two denominators, massless 



Here we present the e-expansion of the massless angular integral with two denominators, 
^j,k{i', 0, 0; e), for the specific values of j, k = —2, —1, 1 and 2. More precisely, we extract a factor 
of J dili_2e and define 



Lk{v;e) = 2-'+''7T' 



r(l - 2e) 



(C6) 



r(i - e) " ^" 

This choice of normalization keeps the expanded expressions simpler and allows for a straight- 
forward comparison with Refs. |3|. Since the results are clearly symmetric in j and k, we can 



24 



restrict to the cases where j > k. We find 
27r 



I-2-2iv;e 
I-i-2{v;e 
h-2{v;e 
h,-2{v;e 
I-i,-i{v;e 
h-iiv;e 
h-i{v;e 
h,iiv;e 
h,i{v;e 



27r 
27r 
27r 
27r 
27r 
27r 

TT 



16 16 80 / 596 656 368 2 \ ^, 

y - + 15" + - 75-"+ 225" 1' + ' 



^ + 1 + 2i; - 6i;2 + 2(^1 + 3i; - 7i;2 



)e + 0(e2) 

^^^^"^''^ + 1 _ 8^; + 6i;2 + 2(1 - 8^; + 9t;2)e + 0(e2) 



4 2 /26 16 . 2^ 



^ + 1 - 2i; + 2(^1 - 2i;)e + 0(6^ 
t; + t;e + 0(e^) , 



e 

1-2?; 
2e 

1 , 1 
— h in — - 



ln2t; + 2Li2(l --y) e + 0{e^] 



TT f 1 „ , 1 

2^ -"' + ^ + ^^^+2 



2i; + 41ni;-ln2?;-2Li2(l-i;) e + 0(e' 



TT 

2^ 



2-f; 



5 + 4t; + (2 - i;) In 



4-6t;-(10-4i;)lni; + (2-t;)ln2?; + 2(2-t;)Li2(l-t;) e + 0(e2)L (C16) 



(C7) 
(C8) 
(C9) 
(CIO) 
(Cll) 
(C12) 
(C13) 
(C14) 
(C15) 



3. Two denominators, one mass 

Next, we present the e-expansion of the angular integral with two denominators and one mass, 
%,fe('yi2,i'ii,0;e), for the specific values j,k = —2,-1,1 and 2. As before, we extract a factor of 
J df2i_2e and set 

■ l+2e_er(l-2e), 



^TT^ _ ^j,k {vi2 , -yii , 0; e) . 



(C17) 
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When A; < 1, the integral is finite in e and we have 



-fi2,-2(^i2,^ii;e) = 27r 



16 8 16 8 2 ^ 

■■"11 - -^Vi2 + ■7TVi2 + 0(ej 



4^-2(^12, ^11 ;e) 



5 5 



2--ui2 + 0(e) 



27r 



15 



(1 - 4v l^^^l ~ ^■"11 ~ lO'^ll + ^'"12 + 16^11^12 - et^is) 

- 2(61;?! - 6i;iit;i2 + vl^ + 2vnvl2) In Q ~ ^| ~ j + 0(e) 



4^-2(^12, ^^ii;e) 



27r 



■uii(l - 4?;ii)V2 



VI - 4i;ii ^^11 + Svji - 12i;iii;i2 + + S^ni^ia) 



+ 2vu (Svu - 2vi2 - 4vuvi2 + SVI2) In Q ^ ^| _ + 0(e) 



/i2-i('"i2>'"ii;e) = 27r 



^i\'_i('yi2,'yii;e) = 27r 



4!-i(^i2,^^ii;e) 



2 - ^ + ^^12) +0(e) 



4 2 
3 ~ 3 

27r 



vi2 + 0(e) 



(1 - 4t;ii)3/2 

vr^(i - 2.12) + (2.11 - .12) In { \~^J,ZZu ) + 



(C18) 
(C19) 



(C20) 



(C21) 
(C22) 

(C23) 



(C24) 



4!-i(^i2,^^ii;e) = - 



TT 



vii{l — 4vii)^/^ 

Vl-4.vn(2vu - vi2) + vu (l - 2vu) In Q~^|~^^|| ^ + 0(e) , 



(C25) 
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On the other hand, for k > 1, the integral has a pole in e. We find 



l'Lli{vi2,vii;e 
7-(i)/ 

-fi2,2(^i2>'yii;e 

-^-i,2(^^i2>^ii;e 
ixX{vv2-,vxi:,e 



I^\{vv2-,vx\:,e 



27r 



27r 



2^12 
TT 



2t; 



12 



+ 1 - 2vxx + 2t;i2 - ^v\2 + 0(e) 



4i;iii;: 



^ + 1 - 2i;i2 + 0(e) 

-l-.(|).o,, 



- 2i;ii + ^12 - uii In ( ^ ) + 0(e) 



12 L 



^12 



= 27r 



27r 



2(.n-2.,2 + 3.f2) + 1 + 4^^^ _ 8^^^ + ^^2^ + o(e)' 



TT 

4^ 



TT 



1 - 2i;i2 
2e 

2vxi - f 12 



- vvi + 0(e) 



+ 2(2i;ii - 2i;i2 + + (2t;ii - t;i2) In (^^j + 0(e) 



(C26) 
(C27) 
(C28) 
(C29) 
(C30) 
(C31) 



(C32) 



12 L 



2fii(3i;ii - 2^12 +^^2) , 1^ 2 i« ,2 , in 2 



+ 2i;ii(3vii - 2^;i2 + vi^) In + 0(e 



'12 



(C33) 



4. Two denominators, two masses 



Here, we present the e-expansion of the angular integral with two denominators and two masses, 
%,fc(^i2 5 ^11) ^22; e), for the specific values j, k = —2, —1, 1 and 2. As before, we extract a factor of 
/ dOi_2e and define 



Ij%vi2,vu,V22;e) = 2 ^+^^7r'^Y^— ^%fc(vi2,^^ii,^^22;e) . 



(C34) 



Clearly this expression is symmetric under the simultaneous exchange of j -H- A; and vu -H- ^22, 
thus we can restrict to the cases where j > k. 



27 



We find 

r(2) 



16 8 16 808 16 



/i2 _2('f^i2> ■"11, ■"22; e) = 27r - ■^■uii - ^vi2 + 7^i;f2 - ii'22 + 7i^^ii'y22 + 0(e) , (C35) 

(2) 4 4 

-fli _2(^i2> ^11,^22; e) = 27r 2 - -V12 - -U22 + 0(e) , (C36) 

4^l2('"i2>'"ii,^22;e) = _ .^5^2 Vl - 4^11 (^1 - lOuii + 2i;i2 + 16^11^12 - 6i;i2 - 2^22 



+ 8t;iiU22^ - 2^6^11 - 6t;iit'i2 + ^12 + 2'Uiifi2 + 2i'iif22 - 8t;iit'22^ 
X In 



1 - VI - 4t;ii 



45^12(^12, ^^ii,^^22;e) 



1 + Vl - 4i;ii 
27r 



0(e) 



(C37) 



uii(l-4i;ii)V2 



Vl - (^vu + 8t;f]^ - 12t;iit;i2 + t'i2 + 8?;iit;f2 



+ 4:VuV22 - 16t;]^it;22 ) + 2f n I Swu - 2f 12 - 4uiiwi2 + 8^12 + ^22 



-f_i_i('Wi2,'yii,^^22;e) = 27r 



4 2 
3 ~ 3 



1^12 + 0(e) 



(C38) 
(C39) 



/Ul(^12,m,^22;6) = 



VI - 4i;ii (^1 - 2i;i2) + (2t;ii - via) In f - 



1 - VI - 4t'i] 



+ VI - 4i;ii 



+ 0(e) 



(C40) 



4^11(^^12, Wii,'y22;e) = - 



TT 



-ffi('yi2,'yii,f22;e) 



vn{l - 4i;ii)3/2 
VI — 4^11 ( 2vii - 'U12 ) +vu{l- 2vi2 ] In 



In 



1 - VI - ^vii 
1 + VI - 4i;ii 



+ 0(e) 



V12 + V ''12 - 4''iif22\ 



^V'wL - 4^11^22 L \'U12 - V^^12 - 4i;il'U22 

4\^('yi2,'yii,'y22;e) ^ 



+ 0(e) 



(C41) 
(C42) 



4^11(^^12 — 4i;iiV22)^/^ 

'■"12 + ^^12 - 4i;iii;22\ 



— ■"11 (^'"12 — 2f 22) In 



y^uf^^^4u^]^(^2i;ii - t;i2^ 
+ 0(e) 



■"12 - Vfi2 - 4?;ii?;22 7 



(C43) 
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41(^^i2,'yii,^^22;e) 



vr 



- 4^11^22 (viivf2 + 8i;iiU22 



8i;ii'U22('yi2 - 4i;iii;22)^/^ 

— VlviiVi2V22 + ^^12 '"22 + 4f llVi2'"22 + ^V\iV22 ~ ^^''^ll''^22) 

— 2?;iii;22 ^3^11^12 — 2i;^2 + t'?2 ~ ^viiV22 + 3vi2t'22 — 4t;ii?;i2V22^ 



X In 



Vl2 + \/^^12 ~ 4t'llt^22 
'"12 - ^/vl2 - 4^11^22 



+ 0(6) 



(C44) 



5. Three denominators, massless 



Finally, we present the e-expansion of the massless angular integral with three denominators, 
^i,A;,z(^^i2j ^^13; f23; e)j for the specific values j,k = 1 and 2. As in the previous examples, we extract 
a factor of J dOi_2e and define 



^j,M(^i2>'yi3>^23;e) = 2 ^''"^^^Y^— ^f^j,fe,;('yi2,'yi3,^^23;e) . 



(C45) 



Clearly this expression is symmetric under the permutations of the indices j, k and I. Hence, we 
can restrict to the cases where j > k > I. We find 

'Vi2 + Vi3 + V23 



h,i,iivi2,vi3,V23;e) 



+ ^^12 - - ^23) lni;i2 



^2,l,l(t'l2,^^13,'y23;e) = - _ 2 2 

8uf2^^13U23 



4i;i2Vi3U23 

+ (vi3 - V12 - V23^ lnt;i3 + (^^23 - Vi2 - Vi3^ ]n.V23 + 0(e) 

""12 + ■"13 + 2t'l2^^23 + 2^13^23 " 2i;i2Vl3'y23 " V. 



TT 



2 

23 



+ - 2^;l2^^l3 + V13 + 4i;i2i^23 + 4t;i3'U23 - 6i;i2'yi3'y23 - '"23 
+ (^12 ^ '^'IS - 2l^l2t'23 - 2f 13^23 + 2vi2t^i3'y23 + vj^^ In t;i2 
+ (yi3 - ■"12 - 2'yi3l'23 - '2'Vi2V23 + 2i;i2'yi3^^23 + ^^23) ^'"12, 

- {^12 + ^13 - 2i;i2'y23 - 2'Ui3'U23 + 2?;i2'yi3^^23 + ^^23) lnt;23 + 0(€ 



(C46) 



(C47) 
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-f2,2,l(^12,«13,«23;e) = - 



+ 3ui2f ?3 - 2wf3 + 6f ^3^23 - 6t;i2Wi3t;23 + 3?;i2'yi3 



+ Gvisvl^ - 6?;i2Vi3'yi3 - ^vfa) ^ + 2^^?2 - 3ui2Wi3 + 8?;i2V?3 - 3?^?3 - 3'i;i2?^23 - 6^12^13^23 

+ 6i;i2Vi3'y23 + 13'yi3'y23 - 16'yi2^'?3^23 + 8'yi2?^i3 + 13?^i3^^23 - 16?^i2'yi3'y23 - 3vi3 

+ (^12 - 3i;i2'yi3 + 2'y?3 - 6i;i3f 23 + 6ui2f i3i;23 - 3t;i2i^i3 - 6t;i3i;23 + 6^121^13^13 + 2v^3^ In vu 

- (^12 - 3i;i2l^?3 + 2t'?3 - 6-Ui3f23 + ^Vi2vl-iV23 + 3'yi2?^i3 + 6i;i3i;|3 - 6?;i2-ui3t'i3 - '^V23^ Invis 

- (vf2 + Svuvl^ - 2vf3 + 6i;i3'i;23 - 6i;i2'yi3i^23 - 3^121^13 - 6i;i3i;|3 + Gvuvnvj^ + 2v^3^ \nv23 



+ 0(6) 



-f2,2,2('yi2,'yi3,'y23;e) = 



(C48) 

(2vf2 - 4'uf2'yi3 - 4'Ui2'y?3 + '^'"13 - 4l^l2^'23 + 6i;f2i'i3i^23 



16uf2i'l3?^23 

- 4i;i3i;23 + 6i;i2'y?3'U23 - 4i;i2i^i3 - 4t;i3t;^3 + 6i;i2Vi3t'i3 + ^V23^ ^ + Avf2 - 12i;f2'yi3 

+ 6Wi2f ^3 - 12tii2f ?3 + 4^13 - 12wf2W23 + 6f i2Ul3t'23 + 19?;i2f 13^^23 + 6f I2f ?3W23 

- 12Vi2'W?3'U23 - 12v?3V23 + 19vi2l^?3'U23 + 6Ui2'ui3 + 6'yi2'yi3'wi3 - 12i;^2 '"13^23 

+ Qvl-iV23 - 12t;i2Vi3i;i3 + 6t;i2Ui3i;23 - 121^121^13 - 12i;i3i^i3 + 19i;i2l'l3l^23 + 41^23 
+ (2'yi2 - 4wiVi3 + 4f i2V?3 - 2t;'^3 - 4t)i2f 23 + 6i;i2Ui3'U23 + 4-i;i3-i;23 

- 6i;i2'y?3l^23 + 4^121^13 + 4'yi3^^23 - 6'i;i2'Ul3'yi3 " '^'"23) ^^^^12 

- (21^12 - 41^12 1^13 + 4i;i2V?3 - 2vf3 - 4vl2V23 + 6^12^131^23 + 41^13 1^23 

- 6i;i2i;i3i^23 - 4i;i2i'i3 - 4i;i3?;|3 + 6i;i2fi3'yi3 + 21^13) In i;i3 

- ^2i;f2 - 4i;?2'f^i3 - 4i;i2'i;?3 + 2i;i3 - 4^12 '^^23 + 6^12 ^13 1^23 - 4i;i3i;23 
+ 6i;i2V?3i^23 + 4'yi2i^i3 + 4'Ui3i;|3 - 6t;i2'yi3i'i3 - 21^23) ^^1^23 + 0(e) 



(C49) 
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